Evaluation of relative entropy of entanglement 
and derivation of optimal Lewenstein-Sanpera 
decomposition of Bell decomposable states via 

convex optimization 

M. A. Jafarizadeh"'^'^ * M.Mirzaee"'^ J M.Rezaee"'^ * 

'^Department of Theoretical Physics and Astrophysics, Tabriz University, Tabriz 51664, Iran. 
^Institute for Studies in Theoretical Physics and Mathematics, Tehran 19395-1795, Iran. 
''Research Institute for Fundamental Sciences, Tabriz 51664, Iran. 

February 1, 2008 



* E-mail: j af arizadeh@t abrizu . ac . ir 
t E-mail: mirzaee@t abrizu . ac . ir 
■l-E-maiLkaramatyQtabrizu.ac.ir 



1 



Evaluation convex optimization and relative entropy of entanglement 



2 



Abstract 

We provide an analytical expression for optimal Lewenstein-Sanpera decomposition 
of Bell decomposable states by using semi-definite programming. Also using the Karush- 
Kuhn- Tucker optimization method, the minimum relative entropy of entanglement of Bell 
decomposable states has been evaluated and it is shown that the same separable Bell de- 
composable state lying at the boundary of convex set of separable Bell decomposable 
states, optimizes both Lewenstein-Sanpera decomposition and relative entropy of entan- 
glement. Keywords: Minimum relative entropy of entanglement, Semi-definite 
programming, Convex optimization, Lewenstein-Sanpera decomposition. Bell 
decomposable states. 
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Quantum entanglement is a powerful property that has attracted much attention, since it pro- 
vides new means of communication, such as quantum cryptography [1], quantum teleportation[2] 
superdense coding, and quantum computation. Therefore, the measure of entanglement of 
composite systems becomes crucial [3, 4, 5, 6, 7]. Even though, for pure bipartite states the 
Von-Neumann entropy is a good measure and it is easily calculated [3]; however, to quantify the 
entanglement of mixed states is a harder task. Some measures have been proposed, among them 
one based on relative entropy[8]. This measure needs a minimization procedure[6, 8, 9, 10]. 
In order to obtain the relative entropy of entanglement of a density matrix , the crucial point 
is to find a separable state ps such that minimizes the relative entropy of entanglement. In 
general case, to find Ps is a laborious work. 

Additionally, there is another useful tool to study entanglement, i.e., Lewenstein-Sanpera de- 
composition (LSD) [11]. Among them is the unique optimal decomposition, the one with the 
largest A in LSD. 

Prom the above considerations, we see that there are some specific separable states playing 
an important role in determining the entanglement measure. Then it is meaningful to deter- 
mine the properties of these separable states. 

On the other hand, over the past years, semidefinite programming (SDP) has been rec- 
ognized as valuable numerical tools for control system analysis and design. In (SDP) one 
minimizes a linear function subject to the constraint that an affine combination of symmetric 
matrices is positive semidefinite. SDP, has been studied (under various names) as far back 
as the 1940s. Subsequent research in SDP during the 1990s was driven by apphcations in 
combinatorial optimization [12], communications and signal processing [13, 14, 15], and other 
areas of engineering [16]. Although SDP is designed to be applied in numerical methods it 
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can be used for analytic computations, too. Some authors try to use the SDP to construct 
an exphcit entanglement witness [17, 18, 19]. Kitaev used SDP duahty to prove the impossi- 
bility of quantum coin flipping [20], and Rains gave bounds on distillable entanglement using 
SDP [21]. In the context of quantum computation, Barnum, Saks and Szegedy reformulated 
quantum query complexity in terms of SDP [22]. The problem of finding the optimal measure- 
ment to distinguish between a set of quantum states was first formulated as a SDP in 1972 by 
Holevo, who gave optimality conditions equivalent to the complementary slackness conditions 
[23] . Recently, Eldar, Megretski and Verghese showed that the optimal measurements can be 
found efficiently by solving the dual followed by the use of linear programming [24]. Audenaert 
at al. [25, 26] apphed convex analysis method and Lagrange duality to calculate the relative 
entropy of entanglement for Werner states. Also Lawrence in [27] used SDP to show that the 
standard algorithm implements the optimal set of measurements. All of the above mentioned 
applications indicate that the method of SDP is very useful. 

In this paper, by using the convex optimization method, we find optimal LSD and minimum 
relative entropy of entanglement for BD states. Finally, we show that the separable states that 
minimize relative entropy of entanglement are the same as states which optimize the LSD 
which lie at the boundary of convex and compact set of separable states [28]. 

The paper is organized as follows: 
In section-2 we give a brief review of convex optimization. Section -3 is about the Optimal 
LSD and also by using SDP we find optimal LSD for Bell decomposable (BD) states. In 
section-4 we explain how one can find the minimum relative entropy of entanglement by using 
the Karush-Kuhn-Tucker(KKT) optimization method, we find the minimum relative entropy 
of entanglement of BD states. Finally it is shown that the same separable BD state lying at 
the boundary of convex set of separable BD states, optimizes both LSD and relative entropy 
of entanglement. The paper is ended with a brief conclusion. 
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2.1 Semi- definite programming 

A SDP is a particular type of convex optimization problem [29]. A SDP problem requires 
minimizing a linear function subject to a linear matrix inequality (LMI) constraint [30]: 

minimize V — c^x 

(2-1) 

subject to F{x) > 0, 

where is a given vector, ), and F{x) = Fq + J2iXiFi, for some fixed hermitian 

matrices Fj. The inequality sign in F{x) > means that F{x) is positive semidefinite. 

This problem is called the primal problem. Vectors x whose components are the variables 
of the problem and satisfy the constraint F{x) > are called primal feasible points, and if 
they satisfy F{x) > they are called strictly feasible points. The minimal objective value c^x 
is by convention denoted by V* and is called the primal optimal value. 

Due to the convexity of set of feasible points, SDP has a nice duality structure, with, the 
associated dual program being: 

maximize —TrlF^Z] 

Z>0 (2-2) 

Tr[FiZ] = Q. 

Here the variable is the real symmetric (or Hermitean) matrix Z, and the data c, Fi are 
the same as in the primal problem. Correspondingly, matrices Z satisfying the constraints 
are called dual feasible (or strictly dual feasible if Z > 0). The maximal objective value of 
—TtFqZ, i.e., the dual optimal value, is denoted by d*. 

The objective value of a primal(dual) feasible point is an upper (lower) bound on V*{d*. 
The main reason why one is interested in the dual problem is that one can prove that d* <V*, 
and under relatively mild assumptions, we can have V* — d*. if the equality holds, one can 
prove the following optimality condition on x: 
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A primal feasible x and a dual feasible Z are optimal which is denoted by x and Z if and 
only if 

F{x)Z = ZF(x) = 0. (2-3) 

This latter condition is called the complementary slackness condition. 

In one way or another, numerical methods for solving SDP problems always exploit the 
inequality d < d* < V* < V, where d and V are the objective values for any dual feasible 
point and primal feasible point, respectively. The difference 

r*-d*^ c^x + Tr[FoZ] = Tr[F{x)Z] > (2-4) 

is called the duality gap. If the equality holds d* — V* , i.e., the optimal duality gap is zero, 
then we say that strong duality holds. 

2.2 Karush-Kuhn- Tucker Theorem: 

All problems treated in this paper are convex optimization problems: minimizing a convex 
objective function subject to upper bound inequality constraints on other convex functions. 
Lagrange duality theory is also well developed for convex optimization. For example, the 
duality gap is zero under mild technical conditions such as Slater 's condition that requires the 
existence of a strictly feasible point. 

We seek first order conditions for local optimality for the following formulation of nonlinear 
programming: 

minimize f{x) 

subject to h{x) — (2-5) 
9{x) < 0, 

where x is a vector in R" . h maps into R'^, and g maps into i?™. 

Kcirush-Kuhn- Tucker (KKT) Theorem: Assuming that functions Qi, hi are differen- 
tiable and that strong duality holds, there exists vectors ( E R'^, and y e R^, such that the 
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gradient of dual Lagrangian L{x* , (^* ,y*) = f{x*) + J2iCihi{x*) + J2iyigi{x*) over x vanishes 
at X* : 

hi{x*) = (primal feasible) 
gi{x*) < 0(primal feasible) 

y* > 0(dual feasible) (2-6) 
y*9iix*) = 
V/(x*) + ^CV/i,(x*) + ^y;Vg,{x*) = 

i i 

Then x* and {C*,y*) are primal and dual optimal, with zero duality gap. In summary, for 
any convex optimization problem with differentiable objective and constraint functions, any 
points that satisfy the KKT conditions are primal and dual optimal, and have zero duality 
gap. Necessary KKT conditions satisfied by any primal and dual optimal pair and for convex 
problems, KKT conditions are also sufficient. If a convex optimization problem with differen- 
tiable objective and constraint functions satisfies Slater's condition, then the KKT conditions 
provide necessary and sufficient conditions for optimality: Slater's condition implies that the 
optimal duality gap is zero and the dual optimum is attained, so x is optimal if and only if 
there are (C*, y*) that, together with x, satisfy the KKT conditions. 
Slater's condition: Suppose x* solves 

minimize f{x) (2-7) 

gi{x) >bi , i = 1, ...,m 
and the feasible set is non empty. Then there is a non-negative vector ( such that for all x 



L{x, C) = fix) + e{b - g{x)) < f{x*) = L{x\ Q). 
If in addition /(.),gij(.) , i — 1, ...,m are continuously differentiable 



(2-8) 



Evaluation convex optimization and relative entropy of entanglement 
In the spatial case the vector x* is a solution of the linear program 

minimize c^x 



8 



subject to Ax = b (2-10) 

x>0, 

if and only if there exist vectors ( e R'^, and y e i?"* for which the following conditions 
hold for {x,c,y) = {x*,C,y*) 

Ax = b (2-11) 

Xi>0; yi>0; x^yj = , z = l,...,m. 

A solution {x*,C*,y*) of (2-11) is called strictly complementary, if x* + y* > 0, i. e., if there 
exists no index i e {1, m} such that x* = y* = 0. 

3 Lewenstein-Sanpera decomposition 

According to LSD [11], any bipartite density matrix p can be written as 

P^Xps + {l-X)pe, Ae[0,l], (3-1) 

where Ps is a separable density matrix and Pe is an entangled state. The LSD of a given density 
matrix p is not unique and, in general, there is a continuum set of LSD to choose from. The 
decomposition with the largest weight A of the separable part is the optimal LSD [36, 37], 
which can be determined uniquely. 

Furthermore, in the case of a two-qubit state it can be shown [11] that the entangled part 
Pe of the optimal LSD is always a pure state i.e, 

p=(l-A)|V')(^| + Ap: (3-2) 
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3.1 Optimal Lewenstein-Sanpera decomposition for Bell-decomposable 
states via convex optimization 



A BD state is defined by 



1=1 



i=l 



where \ipi) are Bell states given by 

1 



V2 



(|00) + |11)), |V'2)= 



V2 



(|oo)-|ii)), 



IV's) = i^n = ^(|oi) + |io))' 1^4) = r ) = ^(|oi) - |io)). 



V2 



In terms of Pauli's matrices, p can be written as, 



V2' 



where 



1 

p = -(/ ® / + Y^tiai (g) CTj), 
i=i 



^1 =Pi -P2+P3-P4, 

t2 = -Pl +;?2 +P3 -P4, 

^3 =Pi +P2 -pa -P4- 



Prom positivity of p we get 



(3-3) 



(3-4) 



(3-5) 



(3-6) 



(3-7) 



1 + tl - ^2 + ^3 > 0, l-ti+t2+t3>0, 
1 + tl + t2 - ts > 0, l-ti-t2-t3>0. 

These equations form a tetrahedral with its vertices located at (1, —1, 1), (—1, 1, 1), (1, 1, —1), 
(—1, —1, —1) [38]. In fact these vertices denote the Bell states given in (3-4). 

On the other hand p given in (3-5) is separable if and only if tj satisfy Eq. (3-7) and. 



1 + tl + ^2 + ^3 > 0, 
l + ti-t2-t3>0, 



l-tl-t2 + t3>0, 
l-ti+t2-t3>0. 



(3-8) 
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inequalities (3-7) and (3-8) form an octahedral with its vertices located at Of = (±1, 0, 0), 
of — (0, ±1,0) and Of — (0,0,±1). So, tetrahedral is divided into five regions. Central 
regions, defined by octahedral, are separable states {pk < |). There are also four smaller 
equivalent tetrahedral corresponding to entangled states(pjt > | for only one of k = 1,...,4), 
where Pk — \ denote to boundary between separable and entangled region. Each tetrahedral 
takes one Bell state as one of its vertices (see Fig-1). 

Now in order to obtain optimal LSD via convex optimization of entangled BD state given 
in (3-3), with pi > |, we first choose an arbitrary separable state as 

4 ^1 

Ps = Y^Pi \<t>i) {<i>i\ , J2p'i = 1, p'l < 7^ (3-9) 
1=1 i=l ^ 

in the separable region. Then according to strict SDP optimization prescription of section 
(2), we need to optimize Tr(Ap'^) with respect to p — Ap'^ > 0, where the feasible solution 
corresponds to 

Amax = "7,zn{ — , — , — , — j. (3-10) 

Pi Pi P3 Pi 

One can show that the only possible choice of A^ax consistent with positivity of p — A^axP's is 

Amax^ ^^~^l\ and p^ = -^,i = 2,3,4. (3-11) 

(1-Pl) Ar, 



'-max 



The equation (3-11) indicates that A^ax is a monotonic increasing function of p[ and its 
maximum value corresponds to p[ — ^, with 

A^„, = 2(l-pi) and pH on ^ = 2,3,4, (3-12) 

^IJ- Pi) 

where this separable state lies at the boundary of separable region, see Fig-1. Substituting the 
results that obtained in p — A^axP's we get 

P - AmaxP's = (2P1 - l)|0l)(0l|, (3-13) 

which is a pure states in agreement with theorem (2) of Ref.[ll]. Equation (3-13) with A^ax 
given in (3-12) is nothing but the optimal LSD of BD states. 
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4 Evaluation of relative entropy of entanglement 



4.1 Minimizing relative entropy via KKT method 

Fundamental to our understanding of distinguishability is the measure of uncertainty in a given 
probability distribution. The uncertainty in a collection of a probability mass function (pmf) 
u> G TZ^ is given by its entropy, 

5(p) = -X:^aog(a;0 (4-14) 

i 

called the Shannon entropy. This measure is suitable for the states of systems described by the 
laws of classical physics, but it will have to be changed, along with other classical measures, 
when we present the quantum information theory. We ultimately wish to be able to talk 
about storing and processing information. For this we require a means of comparing two 
different probability distributions, which is why we introduce the notion of relative entropy 
(first introduced by KuUback and Leibler, 1951 [32]). 
Definition of the relative entropy: 

Suppose that we consider the problem of estimating a probability mass function (pmf) 
00 G TZ'' given a strictly positive prior q G TZ^. The relative entropy between these two 
distributions is defined as [8] 

k 

X{u;q)^Y.^dog- (4-15) 

1=1 9i 

This function is a measure of the distance between coi and Qi, even though, strictly speaking, 
it is not a mathematical metric since it fails to be symmetric X(ciJi \ \ qj) ^ I{qj \ \ cOi). 

A common restriction on a; is a mean constraint, so that if there are k observations A — 
[ai,...,afe] G 7?.* and a mean h ^W^, then uj must satisfy Auj — b. A standard approach 
[31, 33, 34] is to minimize the relative entropy function (4-15) over the constrained probability 
simplex 

l^u; = 1, uj>0 

(4-16) 

Auj = b, 
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where the first two constraints as due to the fact that oj represents a distribution and the 
symbol 1 denotes a vector of ones. Consider the relative entropy optimization problem 



subject to constraint (4-16). The optimization problem (4-17) is convex; in other words its 
objective function is convex, and the equality constraints are linear. Under a suitable constraint 
qualification, such as Slaters condition, the first-order optimality conditions of (4-17) are in 
fact both necessary and sufficient (see, for example, [30, 35]). 

One simple constraint qualification is Slater's condition, i.e., there exists an u in the relative 
interior of the feasible set, such that 



Let C and y be the Lagrange multipliers associated with the first and second constraints of 
(4-17), respectively. 

Definition: (KKT optimality conditions): A triple {uj*, C*, y*) is a first-order KKT point 
of (4-17) if it satisfies the following conditions: 



minimize I{uj] q) 



(4-17) 



a; > 0, l^o; = 1, Au ^ b. 



(4-18) 



(4-19) 



V,X(a;, g) + CI + A^y = 0. 



Theorem: Suppose that Slaters condition holds. Then the vector with components 




(4-20) 



for 



j — 1, /c, solves (4-17), where 



Uj = qjexp{-ajy*), 



(4-21) 



and y* is the Lagrange multiplier corresponding to the constraint Aco — b [10]. 
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4.2 Evaluation of relative entropy of entanglement of Bell decom- 
posable states via KKT method 

In order to obtain minimal relative entropy of entanglement we choose entangled BD state 
(given in (3-3)), with pi > \ and then we choose an arbitrary separable states 



p's = Ylp'il'^i 



1=1 



4 ^ ,1 

0<Pi<l, Y.Pi = 1, p[ = 2> 



(4-22) 



1=1 



at the boundary of separable region. 

Then according to the prescription of the minimization of the relative entropy given in 
section (2), all we need is to minimize I{uj; q) with u'^ — {p' i,p' 2,p' ^/p' = {pi,P2,P3,P4:) , 
Pi < I and Pi > I subject to the constraints (4-16)with matrix A and vector b of the following 
form: 



'^ 1 ^ 









= ( 6i ) , 



(4-23) 



where bi < |. Now, using the solutions (4-20) and (4-21) of KKT (optimality conditions) or 
relative entropy optimization given in (4-19) we obtain 



=piexp(-y*) , U2^P2 , u-i^ps , u^^p^, 



(4-24) 



now, using the slater's constraint we obtain 

Pi exp(-?/i; 



(4-25) 



j-i — - - 

Pi exp{-yl) +P2+P3+P4 

After some calculation one can find the Lagrange multiplier yl and the other components of 
optimal u!*, where we quote only the optimal values uj* , j = 2, 3, 4 in the following: 



, 00*2 = , ^ = (S)/>3 , u;2 = ClE^)p, . (4-26) 



-pi 



pi' 



-pi- 
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Finally inserting to optimal lu* in (4-26), we obtain the following expression for the relative 
entropy: 

T(a;*; q) = 0;^ log ^ + u; log ^ + a;* log ^ + log ^ = h log(^) + (1 - h) log(i^). 

Pi P2 PS Pi Pi 1 - Pi 

(4-27) 

The equation (4-27) indicates that T{uj] q) is a monotonic increasing function of hi and its 
maximum value corresponds to Pi — bi — ^, with 

x(-;<;) = ilog(^) + ilog(^(^). («8) 

and 

The BD separable state with p'- = u* = 2(i-pi) ^^ich minimizes the relative entropy of 
entanglement of entangled BD state p — Yh=i Pi < \ > with pi > | is the same as the one 
that has been found by Vedral and Plenio in Ref. [5]. Also one should notice that c = 2pi — 1, 
where c is concurrence of BD states and we have 

X{u-q) = -Uog{l-^), (4-30) 

which is in agreement with Ref. [6]. 

Finally comparing the results of sections 3 and 4 we see that the same separable BD state 
lying at the boundary of convex set of separable BD states, optimizes both LSD and relative 
entropy of entanglement. 



5 conclusion 

Here in this work, using the convex optimization method, we have been able to obtain optimal 
LSD of BD states ( by using the semidefinite programming method) and we have also evaluated 
the relative entropy of entanglement of BD states with respect to convex set of BD separable 
states ( by using KKT method). 
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At the end it is shown that the same separable BD state lying at the boundary of convex 
set of separable BD states, optimizes both LSD and relative entropy of entanglement. 
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Figure Captions 

Figure- 1: All BD states are defined as points interior to tetrahedral. Vertices Pi, P2, P3 
and P4 denote projectors corresponding to Bell states given in Eqs. (3-4) . Interior points of 
Octahedral correspond to separable states. 



